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We derive an exact expression for the differential conductance for a quantum dot in an arbitrary 
magnetic field for small bias voltage. The derivation is based on the symmetric Anderson model using 
renormalized perturbation theory and is valid for all values of the on-site interaction U including 
| the Kondo regime. We calculate the critical magnetic field for the splitting of the Kondo resonance 

to be seen in the differential conductivity as function of bias voltage. Our calculations for small 
field show that the peak position of the component resonances in the differential conductance are 
reduced substantially from estimates using the equilibrium Green's function. We conclude that it 
is important to take the voltage dependence of the local retarded Green's function into account in 
interpreting experimental results. 



in 
o 



PACS numbers: 



I. INTRODUCTION 



Measurements of the temperature dependence of the differential conductance dl/dVds through a quantum dot 

g 

o 



in the linear response regime have been shown to be in good in agreement with theoretical predictions based on 
magnetic impurity models Measurements of dl/dVds through a quantum dot in the presence of an applied 

£j ' magnetic field, however, have revealed an apparent discrepancy between the magnetic splitting in the Kondo peaks 
'— 1 1 deduced experimentally and theoretical predictions based on the same models 0, S 13 • The differential conductance 
1 ■ measurements are made in the presence of a finite applied bias voltage Vds, which is a situation for which there are 
few theoretical predictions. There are calculations based on the non-crossing approximation (NCA) 6J. This method, 
T^-j- . however, is known to break down in the Fermi liquid regime in the equilibrium case, so there is some uncertainty 
[~^. ' about the validity of the predictions in this case. There are also some recent perturbation calculations of Fuji and 
7-H , Uedo using the Keldysh formalism, taking diagrams up to fourth order in the on-site interaction U into account. 

This approach should reliably describe the trends in the approach to the strong correlation limit, as the interaction 
strength is increased, but cannot fully describe the Kondo regime, where there is an exponential renormalization of 
the Kondo resonance width as a function of U. Rosch et al. have calculated the differential conductance using a 
perturbative renormalization group calculation for the exchange coupling J to leading order in l/ln(Vd s /Tk), which is 
appropriate for the very high voltage regime (Vd s ^> Ik) (see also H3)- This approach is, however, not suitable 
for investigating the Kondo regime, where Vd s and the magnetic field are of order of the Kondo temperature Tk. 

More accurate methods || |3, El El f° r the Kondo regime exist for calculations in an applied magnetic field, 
which have been used to interpret the experimental results 0. 01 ■ but they are restricted to the linear response 
(equilibrium) regime. It is not clear that using the equilibrium response for dI/dVd s in the presence of an applied 
bias voltage will lead to reliable predictions. In this paper we examine the effect of a finite applied bias voltage in 
the small bias regime. In the limit of a small magnetic field B, we calculate asymptotically exactly the shifts of the 
peaks of the differential conductance of particular spin component for small, but finite Vds- The results are based 
k^J ' on renormalised perturbation theory and numerical renormalization group (NRG) calculations. We also deduce the 
value of the critical field B c for two distinct peaks to be seen in dI/dVd s - The knowledge of dI/dVd s to order V^ s for 
' arbitrary values of B is sufficient to determine B c exactly. 

In our treatment the quantum dot will be described by a single impurity Anderson model with particle-hole sym- 
metry coupled via left and right leads to two reservoirs of free electrons. The Anderson model has the form, 

Ham = £ dadld a + Un d ^n d ,i + X](^4 c fc, CT + V k c l,a d a) + Y e fc c iL c fe,o-= C 1 ) 

where tda = £d — ^9^bB /2 is the energy of the impurity levels with spin a = ±1 in a magnetic field B, U is the 
interaction at the impurity site, and Vu the hybridization matrix element to a band of conduction electrons with 
energy e^. In the wide band limit the hybridization weighted density of states, A(cj) = 7r^ fe |Vfc| 2 <5(a-> — e^), can be 
taken as a constant A. Taking this limit justifies the neglect of any magnetic field term acting on the conduction 
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electrons in (Clogston- Anderson compensation theorem [T^)- Any polarization of the conduction electrons only 
affects the impurity via the hybridization function A(lu) and any change to the conduction band density of states due 
to an applied magnetic field is only at the band edges (±-D), and hence negligible if g/isB <C D. Therefore, if A(lo) 
is independent of lu the magnetic field does not need to be included in the electron band. 

The single impurity model can be applied to a quantum dot which is equally coupled via leads to left and right 
reservoirs. In this situation the charge on the quantum dot is coupled only to an even combination of states from the 
left and right channels by an effective hybridization Vk, so it can be mapped into a single channel model. If there is 
a potential difference eVds due to a bias voltage Vds (e is the electronic charge), then the chemical potentials in the 
left and right reservoirs differ from the average chemical potential \i by eVd s /2 and —eV ds /2, respectively. 

A general expression for the current I due to the applied voltage through a quantum dot has been given by Hershficld 
et al. and Meir and Wingreen 0,0, and specialising it to this symmetric case it takes the form, 

oo 

I=C -^T. J MfLn-f R (u)][-lmG d ^,V ds )], (2) 

a — oo 

where G d (co, Vds) is the steady state retarded Green's function on the dot site, and Jl{oj), are Fermi distribution 

functions for the electrons in the left and right reservoirs, respectively, and h is Planck's constant. 
The retarded Green's function G da .(u>, Vds) can be written in terms of a self-energy ^{ui, Vds)-, 

G d „(iv, V ds ) = ./ v - ( v : (3) 

lu - e d + i A - S ff (w, V da ) 

Oguri p^.ll9| has shown that this Green's function for the symmetric model (e d = —U/2) in the Fermi liquid regime 
in the absence of a magnetic field to order to 2 and (eV ds ) 2 , can be expressed in the form, 

G da {u,V ds )= - ' (4) 
oj + iA ~ h,{u), Vds) 



where S(o>, V^ s ) is a renormalized self-energy |20| given by 

S(w, V ds ) = -ic 



2 + \{eV ds f 



with c = [ | (5) 
2A UA 



expressed in terms of the renormalized parameters A and U . These parameters |20j are defined in terms of the local 
impurity self-energy £(cj) and local irreducible 4- vertex Tn(uii, u>2, U4) 20] for V ds = 0, 

A = zA, L> = z 2 r u (0, 0,0,0), (6) 

where z is given by z = 1/(1 — S'(0)). The renormalized parameters A and U, can be calculated directly from NRG 
calculations, as described in references [22I |23| . where they are given in terms of the 'bare' parameters of the model 
A and U. Alternatively they can also be deduced from the Bethe ansatz results po| . 

We can generalise this result by including a magnetic field term [2H. |24| to lowest order in B, and then G drT (uf) 
takes the form, 

G^M = . (7) 

lu + crfjb + iA- V ds ) 

where a = ±1, b = g/iB-B/2, and 77 is a renormalization parameter defined as rj(b) = z(b)(l — £|(0, b)/b). This 
renormalized parameter can be deduced from NRG calculations as shown in reference |23l |25| . We give an example in 
the Kondo regime for U/irA = 4 in figure ^ where the ratio of the parameters to their values in zero field are plotted 
as a function of log(6/TK). As b — > 0, fj(b) is independent of b such that fj(b) — > R, where R is the Wilson or 'x/7' ratio 
for zero magnetic field. In terms of the renormalized parameters for finite field R(b) is given by R(b) = 1 + U (b)p d (0, b) 
[23,|23|, where the free quasiparticle density of states p da (uj, b) is given by 

p°>,&) = A(b)/ " . , (8) 

' (w + abf,{b)f + A 2 (&) V ^ 
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FIG. 1: b dependence of the renormalized parameters for strong coupling (U/nA = 4), in relation with their b — values, 
U = tv A = 4T K = 0.0817A. Note the logarithmic scale of the y-axis. 



which is independent of the spin index a for w = 0. There are in general B 2 and higher order terms in the self-energy 
but for the moment we work only to first order in B so the renormalized self-energy Y> a (uj,Vds) is still given by 
equation JSJ) with 

which reduces to JSJ for b — > 0. 



II. WEAK FIELD LIMIT 



We look at the effects of the different contributions to dI/dVd s in turn. We consider the simplest situation to begin 
with and ignore the renormalized self-energy term in J7J). The differential conductance J5J) at T = then takes the 
form, 



dl 

dVds 



E 



(eVds/2 - lyfjb) 2 + A 



2 ' 



(10) 



where v = ±1. The total result is expressed as a sum over v = ±1, rather than as a sum over a = ±1, as both 
spin up and spin down states contribute to each resonance. Note the result does not simply correspond to the non- 
interacting case, as it includes the many-body renormalization factor f). It corresponds to the case in which the 
interaction between the renormalized quasiparticles (U) is neglected. In the Kondo regime rj = 2, so the effective 
Zeeman splitting of the resonance in this limit is twice the Zeeman splitting (26) for non-interacting electrons (U = 0). 
Note that the energy level on the dot is at the average chemical potential /i, which is shifted by eVd s /2 with respect 
to the value at zero bias. The individual terms in (|10|) for each v have maxima at ±2rjb. The maxima of the sum in 
general occur at a shifted position eVj: = ±2fjbf c (rjb, A) with 



fc(fjb,A) 



1-1- 



A/fjb 



1/2N 



1/2 



(11) 
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provided that fjb > A/V3- For fjb ^s> A can be approximated by f c (fib,A) ~ 1 — | . The peak splitting is 

therefore e(V d + -V ds )= 4ijbf c (fjb, A). A maximum of the differential conductance, therefore, occurs when one of the 
quasiparticle peaks is coincident with left Fermi level at [i + eVds/2 and at the same time the other peak coincides 
with the right Fermi level, /i — eVds/2. This is illustrated in figure|21 In the interpretation of the experimental results 



H + eV da /2 




M - eVds/2 



(a) 

FIG. 2: A schematic plot of the spectrum on the dot and chemical potentials for left/right lead (fj, ± eV/2) and dot (/j) for (a) 
zero bias and zero magnetic field and (b) finite voltage and finite field. 



of dl/dVds the splitting of the Kondo resonance AKondo was identified with the voltage splitting e(V^ s — V^") = 
AKondo |4j, |5j . However, for the renormalized quasiparticles the Kondo splitting is Axondo = 2ijb, which substituted in 
the above result gives e(V^ s — V^~) = 2AKondo- The factor of 2 difference is due to the many-body shift in the energy 
level on the dot, which was not taken into account in the interpretation of the experimental results in 0,01 . 

We should also take into account the terms in the renormalized self-energy T,(u>,Vd s )- For the moment we ignore 
the Vds dependence of Vds), and then 



dl 

dVds 



e 2 A v 
h ^ 



(uj - vfjb) 2 + (A - S J (w)) 2 



(12) 



sW2 



where & (oj,Vd s ) is the imaginary part of Vds) as given in equation JSJ. If the component spectral densities 
(y ± 1) have maxima at u>± then this expression has maxima as a function of Vds for 2ui±. Working to lowest order 
in 6, we find from the exact result for the self-energy to order w 2 , that for a maximum, 



u> = ±fjb - 



A <9£ 7 

~2~"a7 



= ±fjb - u;Ac{b) 



The peak position in the resonances in the spectral density to first order in b is given by 



±f]b 



1 + Ac(b) 



(13) 



(14) 



For b 



find Ac(b) — > | (^) = |(-^ — -^) 2 anc ^ so tms resu lt corresponds to that of Logan and Dickens @, and as 
derived from the renormalized perturbation calculations [2l| . It shows that the effective Zeeman splitting of the free 
quasiparticles, 2?y6, is reduced in the spectral densities by a factor arising from the (J 2 term in the imaginary part of 
the self-energy. In the Kondo limit, U = ttA = 4Tk, so U/irA = 1 and 77 = R = 2; hence the peak position of 2b is 
reduced by a factor of 2/3. This effect will be reflected in the voltage difference between the differential conductance 
peaks, 



<y d + s -v d -) = 



4rjb 



1 + Ac(b) 



Uvb,A), 



(15) 



where f c is a correction factor due to the overlap of resonances similar as above, which can be computed numerically. 
However, in calculating this voltage difference we cannot neglect the voltage dependence of the self-energy. 
When we take the Vd s dependence of the retarded Green's function into account, 



dl 

dVds 



e 2 A^ 



ImG dtT (eVds/2,V d s) 



2eA 
~h~ 



E 



eV d s/2 



du> 



-Im 



dGd,<r(uJ, V ds ) 



dV„ 



ds 



(16) 




FIG. 3: The shift of the component resonance {v — 1) in the differential conductance (units of e 2 /h) in a magnetic field for 
b/A = 0.05 as a function of the bias voltage eVds/A, according to the inclusion of different contributions as described in the 
text. The arrows indicate the respective maxima. 



In considering the contribution from the first term on the right hand side we can replace the eVd s in the self-energy 
by 2u>, and then see where it has a maximum as a function of to. This will correspond to the calculation we have just 
done except that the oj 2 term in Tii(u>) will have an extra factor of 4, as we can see from equation (J5J. We get the 
modified result, 

e «-- F -> = rJb) A( * A) ' (17) 

In the Kondo regime and for b — > the reduction factor arising from the self-energy including the voltage dependence 
is now 1/3, so the splitting on including this term is reduced by an extra factor of 2 when the voltage dependence in 
the retarded Green's function is taken into account. We do need to take account, however, of the second term on the 
right hand side of equation (|16|l involving an integral over u>. This term can be written as 

Hb)e 2 v AV r W2 (w + vflb? — (A + c(b)(u, 2 + 3(e^ s ) 2 /4)) 2 
h dS ^Jo ^ [(^ + ^6) 2 + (A + c(6)( W 2 + 3( e y ds ) 2 /4))2]2- 

If we plot all the contributions to dl / dVds in the very weak field regime then, due to overlap, no magnetic field 
splitting can be observed. We can calculate, however, the shifts in the component resonance for v = ±1. In figurcElwe 
plot the terms in the differential conductance (in units of e 2 /h) given by equation (|16(l as a function of eVds/A, where 
we use v = 1 in 1|12|) for the first term in (|16[) and (|18|l for the second term in We take values corresponding 

to the Kondo regime, with R = fj = 2, b/A = 0.05 (irA = 4Tk). We have distinguished between the different 
contributions, (a) is the case for the non-interacting quasiparticles, (b) includes the to 2 term in the renormalized 
self-energy, (c) includes the first term in equation and (d) takes into account the full expression including the 
integral term in JTSJ. We see that the integral term arising from the voltage dependence of Gd, a (u, Vd s ) causes a 
significant further reduction of the magnetic shift beyond that estimated from the first term in equation l|16() . In an 
experimental conductance measurement, this is not observable, however, due to the overlap of the two components. 



III. CRITICAL FIELD 



Our derivation is restricted to the regime where eVds is small compared to A. These results are sufficient for us to 
deduce the critical value of the magnetic field b c at which two distinct peaks begin to appear in the total differential 
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FIG. 4: The total differential conductance (in units of e 2 /h) in the Kondo regime for larger magnetic field values, calculated 
using equation l|16[l taking into account the full selfenergy expansion from 11211 . These results are asymptotically exact for 
eVds/A <C 1 and approximate, based on a second order expansion in eVds for larger values. 



response. For values of b < b c the differential conductance will have a maximum at Vds — 0, and for b > b c this will 
become a minimum. From the coefficient of dl /dVds to order V$ a we can determine the point at which it changes sign 
as a function of 6, and hence determine b c . The contribution from the integral term in equation l|18|) to order V% s can 
be evaluated trivially, as it is sufficient to this order to put lu = Vds = in the integrand. As a first estimate using 
the above results the value of b c can be calculated analytically, and the result expressed entirely in terms of A and 
the Wilson ratio R = 77(0) = 1 + U/nA, 



bl _ V9 + 20(fl-l)2(l + 5(fl-in-3 

A 2 10R 2 (R-1) 2 ' 1 ' 

In the non- interacting case, R = 1 and b c j A = = 0.577, which corresponds to f c (fjb, A) = from equation 

and in the Kondo regime, R = 2, A = ATk/tt, and 6 c /7k = 0.582, with 7k given by 

Tk = A (J£j V2 e -^/8A+-A/2a (20) 

If the Vds dependence of the Green's function is neglected the result in the Kondo regime is & C /Tk = 0.491; significantly 
smaller than if this term is included. 

The estimated critical magnetic field is comparable with A, and for U ^ it may not be sufficient to work to linear 
order in B. It is possible to work with an arbitrary magnetic field, but in this case the renormalized parameters 
become field dependent. The renormalized self-energy to order to 2 and V£ s can be expressed in the form, 



V ds ) = -c{b) 




a w (&)w 2 + ay (6) 




(21) 



where c(b) is given in equation @. The coefficients for the expansion of the real part of E CT (o;, Vd s ) are 



(22) 
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where £(6) = np d (0,b)ij(b)b, 1(b) is the integral 

OO OO 

1(b)- J /g;K)G°K + c')[G?( W ')] 3 ^^, (23) 

— OO — OO 

where [G° (w)] -1 = lo + afjb + iAsgn(w) is the renormalised free propagator for T = 0. a u (b) is obtained from the 
second order derivative of the renormalized self-energy evaluated at u> = using the renormalized perturbation theory 
to order U 2 . This reduces to the exact result of Yamada (2*3 f° r the symmetric model for 6 = 0. The corresponding 
coefficient ocy(b) for finite voltage is calculated from the U 2 contribution to the retarded self-energy in the Keldysh 
formalism [26|. where the propagator Gq ~(lo, Vds) taken to order V d 2 for T = is 

(uj — arjby + l\ z 

where 6'(ui) is the derivative of the delta- function, together with a similar equation for Gj~(w, Vds)- The result for 
a(b) is 



av(6) = l+ A(6) 



2^(^)^(6)6 



— tan 

A(6) 




(25) 



The contribution to the real part of the self-energy to order cu 2 and V 2 S is proportional to crfj(b)b. It might be thought 
that such a term should cancel out in taking the sum over the two spin components. However, there is a er-independent 
contribution from a cross term with the effective Zeeman term abrj(b), which has to be included. In the limit 6^0 
equation i(?T|) reduces to (JSJ. 

The equation for the critical field b c becomes 



b 2 _ y/(3 - ajbMb)) 2 + 4 7 (6)(5 - a(b))(l + 5 7 (6)) - 3 + «(6) 7 (6) 
A 2 (6) 2 7 (6)(5-a(6))f?(6)2 1 } 

where a(b) = a u (b) + ay (6), and 

7(6) = ^A(&)l> 2 (6M(0,6)] 3 = ^(6)^(0,6)^(6) - l) 2 . (27) 

Equation (|26|l is an implicit equation for b c which can be solved by iteration starting from the much simpler result 
(|19fl . obtained within the linear approximation. 

For a strong coupling situation (U/irA — 4) the result for the critical field obtained by iterating equation (|26|l and 
using the 6-dependent renormalized parameters taken from figure 2] is b c ~ 0.459 A = 0.584Xk- This differs only by 
0.3% from the value obtained from l|19l) . The small difference is due to the fact that the various correction terms due 
to the 6 dependence of the parameters in the more general formula l|26l) tend to cancel giving only a small resultant 
change. 

Plots of the total differential conductance for various fields above and below the critical field are displayed in figure 
01 We have taken the full self-energy as given in (|2f (I into account, including the field dependence of the renormalized 
parameters. Our results are asymptotically exact only for small Vds and a more complete theory is required to 
calculate the splitting at larger bias voltages. The major problem to be solved is the dependence of the self-energy on 
the voltage bias term, when eVds is comparable with and greater than the Kondo temperature Ik, so that a detailed 
comparison with experiment can be made with the experimental results in this regime. 



IV. CONCLUSIONS 



Measurements of the differential conductance in a magnetic field have been used to infer the magnitude of the mag- 
netic splitting of a Kondo resonance The results reveal an apparent disagreement with theoretical predictions, 
but the comparisons with theory have been based on calculations of the equilibrium Green's function for the dot. It 
is not clear that this will constitute a reasonable approximation at the finite bias voltages used in the experiment. 
To examine this question we have derived an expression for the differential conductivity through a quantum dot, 
described by a symmetric Anderson model, for small bias voltages and arbitrary magnetic field. This has enabled us 
to estimate the effect of the voltage dependent terms from the non-equilibrium dot Green's function. Our estimates 
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of the shifts of the component resonances for small magnetic field values differ significantly from the values obtained 
using the equilibrium Green's function. We conclude that it is important to take this voltage dependence into account 
and use the non-equilibrium Green's function for a meaningful comparison with experiment. 

Though our calculations are restricted to the small voltage regime we have also been able to estimate the value of 
the critical field B c for the emergence of two distinct peaks in the total differential conductance. Our approach is 
applicable in all ranges of the interaction U from the weak coupling to the strong coupling Kondo regime for T = 0. 
They complement the perturbation calculations in U of Fuji and Ueda for the same model. It is difficult to make a 
comparison with their results, however, because the interaction simultaneously modifies the renormalized parameters, 
such as 77(6) , and introduces a dependence on Vd s , both directly and through the lo dependence in the self-energy. 
These three effects of the interaction can be clearly distinguished within our formulation and taken separately into 
account. It should be possible to extend our calculations to larger values of V^ s using a renormalized version of the 
perturbation theory, and this approach is currently being investigated. 
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